STIC  FIlE  copy 


Copy  12  of  35  copies 


I  AD-A182  501 


IDA  MEMORANDUM  REPORT  M-199 


GLOBAL  EFFECTS  SIMULATION  STUDIES 


Ernest  Bauer 
Frank  A.  Albini 
Craig  Chandler 


DTIC 

ELECTE 

JUN  1  1  1987 


r  i 


D 


V 


November  1986 


DlSTglBUTlQN  STATEMENT  A 

Approved  for  public  releasof 
Diatribution  Unlimitad  ».• 

Prepared  for 
Defense  Nuclear  Agency 


!• 


!W 

a 


INSTITUTE  FOR  DEFENSE  ANALYSES 

1801  N.  Beauregard  Street,  Alexandria,  Virginia  22311 


IDA  Log  No.  HQ  86-31157 


DEFINITIONS 


to  report  tko 


Reports 

Rapcrtt  an  tka  mu  aalkarttattva  and  Mat  carvMhf  CMaMmi  praPacta  IDA  pMhkii. 
Tbay  aarmally  mMy  raaatta  at  najif  prajacts  vMd  (a)  kava  a  Ptract  kaariap  aa  tacWan 
pHocUoi  mojor  pmnmo,  or  (k)  i Mm*  Ittoot  m  rtpfflcool  trorore  to  tko  Encottw 
Braadi,  Dm  Caafcata  ad/ar  Dm  paklte.  ar  (c)  adnat  tsaaaa  tkal  kava  atfaHIcaat  tcaaaaDc 
kapllcatiaaa.  IDA  Raparta  art  rtviawaP  by  aataiPa  paaals  at  axparta  la  aaaara  Ikak  Mfk 
^aaray  m  nmnsaaa  st  in  pfoweRi  oowfii  inn  way  mMiNi  ay  §*»  riviiNaa  a  idr. 


rvjMvv  aanmiy  aaarava  mauvviy  nsncai  an 
tka  rtapita  at  apacial  aaatysaa,  Maria  raparta  ar 
vat.  Papars  ara  rarlavat  la  aaaara  Mat  (bay  mi 
rataraai  papara  la  prvtaaalaaal  laaraala. 


at  a  tori,  ap  kac  ar  patafc  raacttaa 
iarPa  afcailar  la  tkaaa  cxpactaP  at 


Memoraedem  Reports 


DM  Mooronokooi  Roporto  on  moS  tor  Sn  e 
route  i^otioffoi  wort  ton  to  plot  reocftoo 

praap  aai  paaal  adkWati  la  tarwarP  MParani 


MOO  or  tko  tpooooro  or  tkO  OOOfyOtl  to 

Mo  rooalti  of  oootrroi  or  ol  worirtop 
t  lo  oooooeolly  matlyzoS  rrt  ooovol- 


Sm  co oho  of  aa 


Raotaw  of 


Hipirti  to  taMoO  to  (Mr 


Tko  raaatta  of  l^kA  oort  oro  oloo  troforrf  tt  trtrkur  rof  lofar^aaf 


Tkfo  MoawtOMk 
tor  Sm  ooo  mO 


rvoofoooo  of  MorttM  parttat.  Tko 
re  aoO  aaatnaO.  oar  ookfoctoe  to  Hk 


UNCLASSIFIED 


mm mm 


It  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2t  SECURITY  CLASSIFICATION  AUTHORITY 

NA 


El 3M »i i  i  [  j  r«y  \  4 m:  n  .mvj  : m  j  i »i ? :  r » m  rn  u  i  »i*  u  i 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  HARKINGS 


1.  DISTRIBUTION/AVAILABILITY  OF  REPORT 


Approved  for  public  release;  distribution  unlimited. 


4.  PERFORMING  ORGANIZATION  REPORT  NUMSER<S) 

IDA  Memorandum  Report  M- 199 


«4.  NAME  OF  PERFORMING  ORGANIZATION 

Institute  for  Defense  Analyses 


Sc.  ADORESS  (City,  Stale,  and  Zip  Cade) 

1801  N.  Beauregard  Street 

Alexandria,  VA  22311 

8a.  NAME  OF  FUNDING/SPONSOfllMQ 

ORGANIZATION 

8b.  OFFICE  SYMBOL 
(If  ftppltcftbto) 

Defense  Nuclear  Agency 

Sc.  ADDRESS  (CRy,  Stale,  and  Zip  Cads) 

Washington,  DC  20305-1000 

It.  TITLE  (lacluSe  Security  CIseeMcallen) 

Global  Effects  Simulation  Studies 


12.  PERSONAL  AUTHOR(S). 

Ernest  Bauer,  Frank  A.  Albini,  Craig  Chandler 


13.  TYPE  OF  REPORT  13b.  TIME  COVERED 

Final  from  _4zaa_  TO  n/86 


14.  SUPPLEMENTARY  NOTATION 


S.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


7a.  NAME  OF  MONITORING  ORGANIZATION 

DoD-IDA  Management  Office,  OUSDRE 


7b.  ADDRESS  (CITY,  STATE,  AND  ZIP  CODE) 

1801  N.  Beauregard  Street 
Alexandria,  VA  2231 1 


(.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

MDA  903  84  C  0031 


to.  SOURCE  OF  FUNDING  NUMBERS 


TASK  NO. 


PROJECT 

NO. 


T-U2-237 


WORK  UNIT 
ACCESSION  NO. 


14.  DATE  OF  REPORT  (Yaw, 

November  1986 


IS.  PAGE  COUNT 

71 


17. 

COSAT!  COOES  j 

FIELD 

GROUP 

SUB-GROUP  1 

IS.  SUBJECT  TERMS  (Cenllnue  an  reverse  N  rmrmtmmy  and  Identity  by  black  number) 

Nuclear  winter,  forest  fires,  plume  rise,  atmospheric  aerosols, 
absorption,  scattering. 


19.  ABSTRACT  (Canlbwa  an  ravaraa  >  nesssssry  ana  Identity  by  Mack  numbar) 

Thit  documant  reviews  the  initial  risa  of  a  lira  plume  in  Vie  atmosphere,  including  onsets  of  Sis  condsnsatlon  of  antrainad  atmospheric  moisture  The  condensation, 
expansion.  and  avaporakon  of  Sis  malar  flood  are  examined,  as  are  Via  r  ala  Ova  aflacJa  of  scattering  and  abaorpbon  in  radiatian  transport  througn  the  gray-white 
amoks  piuma  and  a  wood  lb*.  For  appboaSona  to  tna  Nudaar  Winter  problsm.  Vie  toSowmg  points  an  mads  as  mcommendadons  lor  further  work  in  cntical  arses: 

1  The  ckmabc  impact  of  a  smoke  piuma  depends  stgnitcanOy  on  in  height  in  lha  abnoaphsn.  which  is  aftoctod  by  tong-lsrm  atmoapharic  motions  and  m  particular 
by  soiar-induoad  Buoyancy  as  was  as  by  tie  initial  pkane  risa.  This  leng-Mmi  behavior  can  be  studied  by  using  Arete  Haze  data 

2  For  a  grven  optical  Vackneee.  Via  abaorptvo  (black,  sooty)  component  of  smoka  (due  largely  to  oil  and  oil-nlaiad  fuels)  is  mufli  mom  aHsctiva  in  mduang  the 
transmsaion  of  sunlight  and  thus  n  producing  ctmatc  cootng  than  Is  Tie  largely  scattering  (whitafgray)  component  which  predominates  in  most  wood  fires.  The 
smoka  characteristics  tend  to  cnange  with  ime.  and  thu#  tor  ctmatc  application  should  be  investigated  as  isle  at  possible  altar  injection  into  the  atmosphere.  This, 
m  turn,  request  using  smoka  bom  (ha  largest  possible  sources  witch  can  be  studied  tor  the  longest  time*.  Explicitly. 

a  For  iorest  Ursa.  Vis  satsltts  analysis  of  Fraser  st  at.  (1986)  has  bean  used  tor  1-  to  3-day  oW  smoke  totoexng  the  very  large  1982  Yukon  liras:  the  results  need 

tobevakdated. 

a  For  oil  Ibas.  even  vw  largest  sources  an  much  smaSsr  toon  tor  wood  (torasVsgriculturai)  liras;  perhaps  the  best  approach  hem  is  to  follow  the  plume  of  an  oil 
pool  *m,  eg.  at  Sandto  Naaonal  Laboratory,  tor  savers!  hours  by  using  a  sampling  abcmfl 

3  Actual  smoka  plumas  are  always  nonurstorm.  even  at  late  tmss.  This  reduces  the  effective  extinction  and  cooling  due  to  a  given  amount  of  smoke  over  the  case 

of  a  uniform  dwtobutton  Experimental  study  of  large  fin  plumes  to  Isle  (met  wd  provide  ■  meatum  of  the  reduction  in  cooling  over  a  uniform  distnbubon  that  is  to 
be  expected. 


20.  DISTRIdUTIONr  AVAILABILITY  OF  A9STRACT 
□  UNCLASS4FIE  Dr  UNLIMITED  IB  9AME  AS 


22s  NAME  OF  RESPONSIBLE  MOtVIOUAL 

Ernest  Bauer 


DO  FORM  1473.  44  MAR 


□  OTIC  USERS 


21.  ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


22b.  TELEPHONE  (IneluAe  Ares  Cede) 

(703)  578-2873 


22c.  OFFICE  SYMBOL 


•3  A PH  t6NI«n  may  bm  um4  until  •ihtilu*. 
AN  •Um*  tNllmi  arm  it— NN 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 

UNCLASSIFIED 


IDA  MEMORANDUM  REPORT  M-199 


GLOBAL  EFFECTS  SIMULATION  STUDIES 


Ernest  Bauer 
Frank  A.  Albini 
Craig  Chandler 


ID A 


INSTITUTE  FOR  DEFENSE  ANALYSES 

Contract  MDA  903  84  C  0031 
Task  T-U2-237 


PREFACE 


This  is  the  final  part  of  work  performed  by  IDA  for  DNA 
under  Amendment  1  to  Task  Order  T-4-237  (see  Appendix  A).  The 
effort  applies  to  the  "Nuclear  Winter"  hypothesis,  and  addresses 
the  partial  simulation  of  this  scenario  by  very  large  forest 
and  industrial  fires,  and  in  particular,  the  observation  of  the 
smoke  plumes  by  using  weather  satellite  imagery.  The  first  and 
major  part  of  this  effort  was  reported  under  M-116  (December 
1985).  It  consists  of  two  items,  namely,  a  panel  report  on  the 
use  of  existing  data  from  weather  satellites  to  study  large 
(forest  and  industrial)  fires  as  experiments  of  opportunity 
and,  as  an  example  of  what  can  be  done,  a  preliminary  report  on 
what  one  can  learn  from  the  December  1982  fire  of  a  large  crude 
oil  storage  tank  at  Tacoa,  Venezuela.  A  revised  version  was 
presented  at  the  1986  IRIS  Targets,  Backgrounds,  and  Discrimina¬ 
tion  Meeting,  and  will  appear  in  the  Proceedings.  The  summary 
of  M-116  is  included  as  Appendix  B. 
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ABSTRACT 

This  document  reviews  the  initial  rise  of  a  fire  plume  in 
the  atmosphere,  including  effects  of  the  condensation  of 
entrained  atmospheric  moisture.  The  condensation,  expansion, 
and  evaporation  of  the  water  cloud  are  examined,  as  are  the 
relative  effects  of  scattering  and  absorption  in  radiation 
transport  through  the  gray/white  smoke  plume  from  a  wood  fire. 

For  application  to  the  Nuclear  Winter  problem,  the  following 
points  are  made  as  recommendations  for  further  work  in  critical 
areas : 

1.  The  climatic  impact  of  a  smoke  plume  depends  signif¬ 
icantly  on  its  height  in  the  atmosphere,  which  is 
affected  by  long-term  atmospheric  motions  and  in 
particular  by  solar-induced  buoyancy  as  well  as  by  the 
initial  plume  rise.  This  long-term  behavior  can  be 
studied  by  using  Arctic  Haze  data. 

2.  For  a  given  optical  thickness,  the  absorptive  (black, 
sooty)  component  of  smoke  (due  largely  to  oil  and  oil- 
related  fuels)  is  much  more  effective  in  reducing  the 
transmission  of  sunlight  and  thus  in  producing  climatic 
cooling  than  is  the  largely  scattering  (white/gray) 
component  which  predominates  in  most  wood  fires.  The 
smoke  characteristics  tend  to  change  with  time,  and 
thus  for  climatic  application  should  be  investigated 

as  late  as  possible  after  injection  into  the  atmosphere. 
This,  in  turn,  requires  using  smoke  from  the  largest 
possible  sources  which  can  be  studied  for  the  longest 
times.  Explicitly, 


vxi 


a.  For  forest  fires,  the  satellite  analysis  of  Fraser 
et  al.  (1986)  has  been  used  for  1-  to  3-day-old 
smoke  following  the  very  large  1982  Yukon  fires; 
the  results  need  to  be  validated. 

b.  For  oil  fires,  even  the  largest  sources  are  much 
smaller  than  for  wood  (forest/agricultural)  fires; 
perhaps  the  best  approach  here  is  to  follow  the 
plume  of  an  oil  pool  fire,  e.g.,  at  Sandia  Nationa 
Laboratory,  for  several  hours  by  using  a  sampling 
aircraft. 

Actual  smoke  plumes  are  always  nonuniform,  even  at 
late  times.  This  reduces  the  effective  extinction  and 
cooling  due  to  a  given  amount  of  smoke  over  the  case 
of  a  uniform  distribution.  Experimental  study  of 
large  fire  plumes  to  late  times  will  provide  a  measure 
of  the  reduction  in  cooling  over  a  uniform  distribu¬ 
tion  that  is  to  be  expected. 
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SUMMARY 


The  Tacoa,  Venezuela,  oil  fire  discussed  earlier  in  this 
study  (see  Part  II  of  M-116)  produces  a  white  plume  which  is 
shown  in  Fig.  1  (p.  1-2).  This  is  interpreted  as  the  capping 
cumulus  cloud  produced  by  the  condensation  of  water  vapor  from 
fuel  combustion  and  the  entrainment  of  moist  near-surface  air 
which  is  transported  up  to  the  plume  stabilization  height  in 
the  relatively  cold  middle  troposphere.  Such  capping  cumulus 
clouds  are  formed  not  infrequently  from  sufficiently  large 
fires  under  appropriate  meteorological  conditions  (moist  and 
warm  lower  atmosphere,  dry  and  cold  upper  atmosphere).  Some 
examples  are  provided  by  large  forest  fires  and  by  large  agri¬ 
cultural  burns  in  the  tropics  (see  Fig.  2,  p.  1-5,  for  example). 
The  direct  climatic  effect  of  the  capping  cloud  is  not  large 
because  it  evaporates  rapidly — at  most  within  several  hours — but 
because  the  latent  heat  of  condensation  of  atmospheric  water 
vapor  in  the  plume  is  released  during  the  plume  rise,  it  will 
cause  the  plume  to  rise  higher  than  it  would  rise  without  this 
energy  input. 

The  climatic  impact  of  a  smoke  layer  depends  on  its  alti¬ 
tude,  with  a  higher  layer  producing  significantly  more  cooling 
than  a  lower-altitude  layer.  The  altitude  distribution  of 
smoke  is  affected  not  just  by  its  initial  injection  but  also  by 
the  long-term  motion  in  the  atmosphere  including,  in  particu¬ 
lar,  its  solar-induced  buoyancy;  perhaps  the  best  way  to  study 
or  simulate  this  is  by  using  Arctic  Haze  as  a  model,  where  on 
occasion  near-surface  emissions  from  a  metal  smelter  at  Norilsk, 
Siberia,  have  been  tracked  for  approximately  ten  days  as  the 
source  of  a  3-km-high  layer  over  Alaska  (in  the  AGASP  program 
see  Schnell  et  al.,  1984). 


Here  a  preliminary  discussion  is  given  of  the  plume  rise 
for  different  types  of  fires,  including  crude  estimates  of  the 
lifetime  of  the  water  plume,  and  also  of  the  optical  properties 
of  the  residual  smoke  plume.  Especially  for  a  wood  fire,  the 
(gray/white)  smoke  plume  mainly  scatters  sunlight  rather  than 
absorbing  it  like  an  idealized  black  soot  plume.  For  a  given 
dimensionless  optical  depth,  the  transmission  of  sunlight 
through  a  scattering  plume  is  much  greater  than  that  through  an 
absorbing  plume,  and  thus  the  "Nuclear  Winter"  cooling  of  an 
absorbing  smoke  (e.g.,  from  oil)  is  much  greater  than  that  from 
a  scattering  smoke  (e.g.,  from  wood). 

Estimates  of  Lawrence  Livermore  National  Laboratory  (Knox 
et  al.,  1986)  indicate  that  the  fuel  in  a  urban/industrial 
location  is  typically  20  percent  oil  and  oil-based  material 
(plastic,  etc.)  as  against  80  percent  for  wood,  paper,  etc., 
so  that  the  scattering  effects  will  be  important. 

Because  smoke  characteristics  change  with  time,  it  is 
important  to  measure  the  late-time  optical  properties  of  smoke 
from  both  wood  and  oil  fires.  This  requires  the  largest  pos¬ 
sible  fires  whose  plumes  can  be  followed  for  a  relatively  long 
time.  (Wall  effects  make  long-term  laboratory  studies  of  smoke 
properties  infeasible  or  probably  unrepresentative.)  Fraser  et 
al.  (1986)  have  studied  the  smoke  plume  of  a  very  large  series 
of  forest  fires  (2.5  million  acres  in  Yukon,  July  1982)  for  up 
to  three  days,  by  an  as-yet-unvalidated  method  which  uses  satel¬ 
lite  imagery.  This  method  should  be  validated  and  extended. 

Oil  fires  are  typically  much  smaller:  the  very  large  Tacoa  oil 
fire  had  a  smoke  generation  rate  perhaps  1/30  that  of  the  Yukon 
fires  (see  M-116)  and  cannot  be  studied  in  this  way.  To  study 
late-time  smoke  from  an  oil  fire,  it  may  be  appropriate  to  use 
the  Sandia  National  Laboratory's  pool  fires  with  a  sampling 
aircraft  which  can  presumably  follow  the  plume  for  several 


By  studying  a  fire  plume  for  an  extended  time,  one  can 
investigate  its  pacchiness  under  representative  conditions. 

This  will  be  important  for  computing  climatic  cooling  because 
an  individual  smoke  plume  is  never  uniform  (although  the  average 
over  many  different  plumes  will  be  uniform).  If  a  given  amount 
of  smok*  is  distributed  uniformly  in  a  horizontal  plane,  its 
extinction  will  be  a  maximum,  so  that  this  nonuniformity  will 
lead  to  a  possibly  significant  reduction  in  cooling.  A  simple 
model  (given  in  Appendix  D  here)  suggests  that  for  optical 
thickness  3  a  reduction  by  perhaps  a  factor  of  1.5  to  8  is 
possible. 

From  the  standpoint  of  their  potential  impact  on  the 
anticipated  climatic  ("Nuclear  Winter")  cooling,  the  following 
recommendations  for  further  work  are  offered: 

1.  Use  Arctic  Haze  data  to  study  the  motion  of  plumes  of 
pollutants  at  very  late  times  (up  to  ten  days). 

2.  Validate  the  satellite  analysis  of  optical  properties 
of  wood  smoke  plumes  at  very  late  times  (days). 

3.  Study  the  optical  properties  of  smoke  from  oil  fires 
(sooty)  as  late  as  possible  (up  to  several  hours) 

4.  Study  nonuniformities  of  smoke  plumes  in  the  atmo¬ 
sphere  at  late  times. 


1. 


INTRODUCTION 


In  Part  II  of  M-116  (Bauer,  1985;  see  also  Bauer,  1986) 
the  characteristics  of  the  smoke  plume  of  the  oil  fire  of  20  to 
22  December  1982  at  Tacoa,  Venezuela,  were  analyzed  and  shown 
to  be  largely  due  to  water  droplets,  evidently  from  the  conden¬ 
sation  of  moisture  generated  by  combustion  and  also  ambient 
moisture  transported  up  to  the  stabilization  height  of  the 
plume.  Figure  1  shows  visible  satellite  imagery  from  NOAA-6, 
and  one  sees  that  the  smoke  plume  looks  white,  like  a  water  (or 
ice)  cloud.  While  the  fire  was  large  in  terms  of  oil  fires 
(the  largest  reported  between  1970  and  1984  on  a  worldwide 
basis),  its  plume  could  not  be  seen  on  satellite  imagery  for 
more  than  five  hours.1 

Table  1  reviews  the  general  characteristics  of  stabilized 
plumes  due  to  large  oil  fires.  Oil  is  well  known  to  burn  with 
a  smoky  flame,  but  nevertheless  the  stabilized  plume  is  some¬ 
times  black  and  sometimes  (near)  white.  For  comparison.  Table  2 
gives  equivalent  information  on  wood  fires.  Smoke  from  intense 
forest  fires  is  often  dense  black  also  (generally  when  combus¬ 
tion  occurs  under  fuel-rich  conditions),  but  usually  the  plume 
turns  white  at  some  altitude,  due  to  the  condensation  of  atmo¬ 
spheric  moisture.  A  capping  cumulus  cloud  is  sometimes  pro¬ 
duced,  and  it  may  survive  for  some  hours. 

For  the  Global  Effects/Nuclear  Winter  application,  the 
main  contribution  comes  from  urban  rather  than  forest  fires. 


^Figure  1  shows  the  plume  at  7:00  a.m.  LT.  GOES  imagery 
at  12:00  LT  displays  a  number  of  smaller  plumes,  while  the 
NOAA  polar  orbiter  overflight  at  2:42  p.m.  shows  so  much 
cumulus  that  the  fire  plume  can  no  longer  be  detected. 
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TABLE  1. 


CHARACTERISTICS  OF  STABILIZED  PLUMES  FROM  OIL  FIRES 


o  Oil  burns  with  a  sooty  flame,  but  the  stabilized  plume  is 
sometimes  black,  sometimes  white. 

o  Oil  storage  tanks  at  Tacoa,  Venezuela  burned  on  20  December 
1982;  the  plume  was  white. *■ 

o  Fire  at  an  oil  refinery  in  Long  Beach,  California  on  22  May 
1958  produced  a  black  plume.2 

o  Meteotron,  S.  France  (large  array  of  oil  burners)  produced  a 
stabilized  plume  that  was  sometimes  black,  sometimes  white. ^ 

o  1  kg  oil  requires  15  kg  air  for  combustion  and  produces  a 
"relatively  large"  quantity  of  black,  sooty  smoke,  also 
~  1.3  kg  water  of  combustion. 

o  During  its  rise  to  stabilization,  the  plume  from  a  high- 

energy  fire  entrains  a  large  quantity  of  near-surface,  i.e., 
relatively  moist,  air  which  it  takes  to  its  stabilization 
altitude.  As  a  crude  estimate,  we  shall  assume  a  100:1 
dilution  of  the  plume  by  entrainment  of  air. 

o  The  atmosphere  at  Tacoa,  Venezuela  would  be  very  moist 
(~  4.5  g/cm2  precipitable  water)  while  that  at  Long  Beach 
in  May  would  be  fairly  dry  (~  1.5  g/cm2)  (see  Appendix  C). 

•  Thus,  1  kg  fuel  produces  perhaps  100  to  200  g  smoke,  includ¬ 
ing  perhaps  10  to  50  g  soot,  and  also  brings  Z  2000  g  water 
(from  combustion  and  entrainment)  to  stabilization  in  the 
atmosphere . 


InOAA-6  satellite  imagery  (see  Bauer,  1985)  who  points  out 
that  for  mean  particle  size  Z  0.1  urn,  the  mass  of  the 
plume  is  significantly  greater  than  a  plausible  mass  of 
soot,  i.e.,  can  only  be  condensed  water. 

2Davies  ( 1959 ) . 

-*L.  Radke  (University  of  Washington),  private  communication. 

4Fresh  soot  consists  of  very  small  (~  0.01  pm)  hydrophobic 
particles  with  a  C/H  ratio  of  about  12/1  by  mass.  The  parti¬ 
cles  tend  to  agglomerate  as  they  age. 
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TABLE  2.  SOME  CHARACTERISTICS  OF  WOOD  SMOKE 


a  Wood  fires  (forest  or  agriculture)  generally  produce  a 
white  or  gray  smoke  plume. 

a  When  a  wood  fire  burns  fuel-rich,  e.g,  as  it  reaches  a 

large  supply  of  fuel,  the  smoke  may  turn  black  momentarily, 
but  generally  it  soon  turns  white  again. 

o  Flaming  combustion  gives  a  relatively  small  quantity  of 
sooty  smoke  while  smoldering  combustion  produces  a  rela¬ 
tively  large  quantity  of  white  or  gray  smoke. 

o  1  kg  wood  requires  approximately  5  kg  air  for  combustion, 
producing  approximately  0.6  kg  water  of  combustion,  1.5  kg 
CC>2,  significant  quantities  of  CO  and  other  gaseous  mole¬ 
cules  including  hydrocarbons,  plus  100  to  200  g  gray/white 
smoke  (made  up  of  complex  organic  molecules,  mainly  liquid 
droplets,  but  including  perhaps  10  to  50  g  soot. 


An  analysis  by  Lawrence  Livermore  National  Laboratory  ( LLNL ) , 
Knox  et  al.  (19B6),  p.  19,  finds  that  available  urban  fuel  is 
typically  15  to  20  percent  oil  or  oil-based  (plastic,  etc.) 
which  would  be  expected  to  produce  black,  sooty  smoke,  while 
the  remaining  80  to  85  percent  is  wood  or  wood-based  (paper, 
etc.)  which  tends  to  produce  white/gray  smoke  that  scatters 
sunlight  rather  than  absorbing  it.  Even  if  the  smoke  genera¬ 
tion  rate  and  particle  size  distribution  are  the  same,  the 
quantitative  cooling  effectiveness  of  white  and  black  smoke 
will  be  quite  different. 

Here  the  problem  of  plume  rise  and  the  condensation  and 
evaporation  of  moisture  in  a  fire  plume  is  reviewed  briefly, 
largely  for  the  case  of  wood  fires.  Figure  2  shows  agricul¬ 
tural  burning  in  Brazil,  where  the  large  fire  (A)  appears  to 
produce  a  capping  cloud,  while  smaller  fires  such  as  (B)  do 
not.  (The  capping  cloud  of  fire  (A)  has  horizontal  dimensions 
of  20  to  30  km. ) 
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FIGURE  2.  Development  of  smoke  plume  from  small  and  large 
fires.  (Sketched  from  space  shuttle  photograph 
I41-D-31-108  of  land  clearing  in  Brazil  in  August 
1984  .  ) 


2.  PHENOMENOLOGY 


2.1  GENERAL  COMMENTS 


Reference  to  Tables  1  and  2  shows  that  the  combustion  of 
1  kg  fuel  yields  0.1  to  0.2  kg  smoke  in  addition  to  -  1  kg  of 
water  of  combustion.  In  addition,  some  1  to  10  kg  atmospheric 
water  vapor  is  entrained  with  the  relatively  moist  near-surface 
air  and  is  taken  to  stabilization  in  the  plume.  Thus,  the  smoke 
plume  always  contains  much  more  water  vapor  than  "smoke".  Dur¬ 
ing  the  plume  rise  to  stabilization  into  a  colder  atmosphere, 
the  excess  water  vapor  may  condense  to  form  a  capping  cumulus 
cloud,  provided  the  fire  is  sufficiently  energetic  and  the  atmo¬ 
sphere  is  sufficiently  moist  (see  Appendix  C  for  a  discussion 
of  atmospheric  moisture). 

The  details  of  plume  rise,  entrainment,  condensation,  and 
stabilization  are  complex  and  are  reviewed  in  Section  2.2  below. 
The  physical  details  of  how  the  condensation  occurs  are  not 
well  understood;  it  probably  will  not  occur  on  soot  particles 
nor  on  large  organic  molecules,  all  of  which  tend  to  be  hydro- 
phobic.  Rather,  it  probably  occurs  on  inorganic  ash  particles 
or  on  sulfuric  acid  droplets  {which  result  from  the  oxidation 
of  sulfur  present  to  a  few  percent  in  oil  or  coal,  but  not  in 
wood).  In  any  case,  there  are  usually  more  than  sufficient 
condensation  nuclei  for  aerosols  to  form. 

The  formation,  spreading,  and  (relatively  rapid)  disappear¬ 
ance  of  a  water  cloud  is  discussed  briefly  in  Section  2.3,  and 
the  late-time  behavior  of  smoke  is  treated  in  Section  2.4. 
Section  2.5  reviews  the  optical  properties  of  smoke  from  both 
wood  and  oil  fires,  raising  some  major  unresolved  questions  on 
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the  composition  of  smoke  and  thus  on  the  relative  contribution 
of  scattering  and  absorption  to  extinction  in  the  visible. 
Finally,  the  climatic  effects  of  smoke  are  briefly  summarized 
in  Section  2.6. 

2.2  PLUME  RISE 

A  large  fire  provides  a  substantial  source  of  energy  with 
mass  injection  (including  water  and  smoke)  on  an  appropriately 
large  geometrical  scale,  which  gives  rise  to  a  (turbulent) 
buoyant  plume  that  entrains  ambient  air  as  it  rises  to  its 
stabilization  height  in  the  atmosphere.  These  qualitative  char¬ 
acteristics  of  the  overall  plume  behavior  are  understood,  but 
as  yet  no  simple,  general,  and  quantitative  description  of  the 
motion  has  been  given.  Thus  Briggs  (1969),  in  an  evaluation  of 
the  rise  and  stabilization  of  chimney  plumes,  reviews  some 
thirty  different  formulas  that  predict  how  effluents  from 
nuclear  power  plants  are  dispersed  into  the  atmosphere.  He 
finds  no  single  simple  but  generally  useful  model  for  the 
process . 

The  flow  field  induced  by  a  source  of  buoyancy,  with  or 
without  mass  injection,  has  been  modeled  by  many  investigators. 
Morton,  Taylor,  and  Turner  (1956)  treated  steady  point  sources 
and  instantaneous  releases  (thermals)  in  uniform  and  uniformly 
stratified  atmospheres.  Their  result  for  the  maximum  height  of 
the  plume,  hmax,  from  a  steady  source  of  heat,  0,  in  a  uniformly 
stratified,  quiescent  atmosphere  is  a  simple  power  law.  For  the 
case  of  a  surface  temperature  of  288  K  and  a  lapse  rate  (rite  if 
decrease  of  temperature  with  height)  6.5  C/km,  they  obtain  the 
resu  1 1 

-  4<>  o‘/4  .  Ill 

where  hmax  is  measured  in  meters  and  0  in  kilowatts. 


Morton  (1957)  extended  this  formulation  to  include  the 
effect  of  moisture  from  the  heat  source  and  also  entrained  from 


the  ambient  environment.  In  this  case,  the  power  law  similar¬ 
ity  solution  leading  to  Eq.  (1)  breaks  down  at  the  height  at 
which  moisture  condensation  occurs.  As  Morton  (1957)  points 
out,  the  fate  of  the  plume  above  the  condensation  height  may  be 
quite  independent  of  the  source,  while  below  it,  the  solution 
for  a  dry  atmosphere  is  an  excellent  approximation.  Thus,  if 
the  plume  reaches  the  condensation  altitude  before  it  comes  to 
rest,  it  may  rise  very  little  further  or,  on  the  other  hand, 
it  may  rise  substantially  higher.  In  many  cases,  its  fate 
above  the  condensation  height  is  determined  more  by  the  struc¬ 
ture  of  the  atmosphere  than  by  the  source  of  buoyancy. 

Near  the  surface  of  the  earth,  air  contains  between  a  few 
grams  and  several  tens  of  grams  of  water  per  kilogram  of  air, 
depending  on  location,  season,  and  weather  (see,  e.g.,  Appen¬ 
dix  C,  Fig.  C.l).  As  the  plume  rises,  it  entrains  and  mixes 
ambient  air  into  the  flow,  slowing  its  rise  and  cooling.  If 
the  source  of  buoyancy  is  dry  (e.g.,  a  source  of  heat  only) 
then  the  moisture  content  of  the  plume  flow  at  each  altitude  is 
due  solely  to  moisture  entrained  at  lower  altitudes.  However, 
the  combustion  of  oil  or  wood  produces  water  approximately 
equal  in  mass  to  the  mass  of  fuel  burned,  while  if  live  vegeta¬ 
tion  is  included,  the  water  due  to  its  combustion  is  about 
twice  that  for  dry  wood.  Thus,  a  fire  plume  will  have  an 
initial  moisture  content  substantially  higher  than  does  surface 
air;  it  is  only  because  it  is  also  hotter  that  this  moisture 
does  not  condense  immediately. 

As  the  plume  rises  into  the  atmosphere,  the  plume  gases 
expand  and  cool.  Under  the  assumption  of  adiabatic  condensa¬ 
tion,  the  Clapeyron  equation  gives  a  relationship  between 
temperature  and  saturation  vapor  pressure  of  water.  This  rela¬ 
tionship  allows  one  to  close  the  set  of  equations  describing 
the  plume  flow,  which  now  contains  a  variable  that  describes 


what  fraction  of  water  is  in  the  liquid  phase.  The  difference 
in  density  between  the  plume  flow  and  the  atmosphere  depends 
on  this  apportioning.  The  density  defect  vanishes  when  enough 
water  condenses  and  the  temperature  in  the  plume  is  low  enough — 
often  several  degrees  below  the  local  ambient.  The  gases  in 
the  plume  follow  a  "wet  adiabat"  on  the  temperature-vs-pressure 
plane,  while  the  moisture  fraction  evolves  according  to  the 
water  entrainment  rate  and  plume  size.  While  the  density 
defect  may  be  very  small  (often  in  parts  per  thousand)  the 
plume  may  be  very  large,  so  that  the  forcing  term  in  the  equa¬ 
tion  of  vertical  motion  can  still  be  quite  large,  as  it  is 
proportional  to  the  horizontal  section  area  of  the  plume  multi¬ 
plied  by  the  density  defect.  Thus,  the  plume  may  rise  very 
slowly,  but  continue  upward  for  a  great  distance  under  the 
influence  of  these  small  forces. 

However,  this  is  not  the  way  the  process  proceeds  in 
nature,  since  the  atmosphere  is  seldom  quiescent.  There  is  a 
wind,  which  generally  increases  in  magnitude  and  changes  direc¬ 
tion  with  increasing  height.  The  change  in  direction  is  usually 
due  to  a  balance  between  shear  stress  and  Coriolis  acceleration, 
resulting  in  a  clockwise  rotation  with  decreasing  altitude  in 
the  northern  hemisphere.  At  midlatitudes  this  deviation  is 
about  45  degrees.  The  wind  field  in  the  near-surface  atmosphere 
is  quite  variable,  and  local  meteorological  soundings  are  neces¬ 
sary  for  accurate  predictions  of  plume  behavior. 

In  fact,  wind  is  usually  the  dominant  factor  in  the 
behavior  of  a  fire  plume.  Indeed,  several  investigators  have 
treated  the  wind-blown  plume  as  embedded  in  the  wind  field  and 
evolving  under  the  perturbing  influence  of  buoyancy.  A  large 
body  of  literature  exists  (see,  e.g.,  Turner,  1960,  or  Briggs, 
1969  for  surveys)  but  the  problem  is  too  complex  to  admit  even 
a  crude  closed  form  approximation.  Briggs  (1969)  discusses  and 
compares  tens  of  empirical  formulas,  finding  none  of  them  con¬ 
vincingly  robust. 
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An  empirical  method  used  to  estimate  the  height  of  plumes 
from  wildland  fires  is  not  discussed  in  either  of  the  sources 
mentioned  above.  This  formula,  given  in  Chandler  et  al.  (1983) 
defines  a  quantity  ATp  in  terms  of  a  mean  power  density 
q  (kW/m2),  the  radius  of  a  circle  that  would  enclose  the  burned 
area,  rQ  (m),  and  the  radius  r^  (m)  that  a  plume  would  have  at 
altitude  h  if  it  grew  linearly  with  height  at  an  angle  a  esti- 


mated 

by  the 

user,  where 
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low-intensity  fires  a  ■■ 
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(2) 

and 

rl 

“  ro  +  h 

tana  . 

(3) 

The  procedure  of  Chandler  et  al.  is  the  following: 

1.  plot  the  dry  adiabat  from  the  surface  to  altitude 

2.  subtract  the  dry  adiabat  from  the  sounding  temperature 
at  altitude  to  obtain  AT(K) 

3.  add  the  ambient  wind  speed  u  (m/s)  to  AT  to  obtain 
(u  +  AT)  ( K) . 

4.  the  plume  tops  out  where 

ATp  =  AT  +  u  .  (4) 

Table  3  lists  some  representative  input  quantities  for 
wildland  fires,  and  also  for  comparison  numbers  for  the  Tacoa, 
Venezuela  oil  fire  and  for  the  "Meteotron"  controlled  oil  fire, 
(Church,  Snow,  and  Dessens,  1980).  Figure  3  shows  some  typical 
profiles  for  AT  and  AT  +  u  gives  the  predicted  height  of 
plume  rise.  In  Figure  3  we  show  Tp  for  the  fire  C  (intense 
brush  or  crown  fire)  from  Table  3.  For  this  fire  the  Morton- 
lay  lor-Turne  r  formula  (1)  gives  a  plume  rise  height  of  5.1  km, 
while  the  Chandler  et  al.  procedure  of  Eq.  (4)  gives  a  plume 


height  of  6.5  km  for  Chandler's  canonical  profile  or  7.7  km 
with  no  wind;  for  the  California  basin  fire  profile,  the  plume 
rise  would  be  5.5  km,  or  for  the  South-East  Asian  Sherwood 
fire  it  is  6.3  km. 

These  numbers  for  Fire  C  are  in  general  agreement  with  one 
another,  and  with  the  predictions  for  the  Morton-Taylor-Turner 
formula  (1)  if  one  assumes  that  the  entire  burned  area  is  alight 
and  burning  with  the  power  density  used  for  the  empirical  calcu¬ 
lations.  Comparable  agreement  is  obtained  for  the  other  wild¬ 
land  fires  A-E.  The  fact  that  these  results  agree  might  give  a 
misleading  impression  of  the  physics.  Recall  that  the  Morton- 
Taylor-Turner  formula  applies  to  a  point  source;  to  apply  this 
to  the  4500  ha  Tillamook  burn  is  unreasonable.* 

To  emphasize  this  point,  the  fires  listed  in  Table  3  were 
represented  by  the  parameters  listed  in  Table  4  which  estimate 
the  actual  area  burning  at  one  time  and  the  power  generated, 
and  computations  were  performed  by  a  numerical  integration  of 
equations  equivalent  to  those  given  by  Morton  (1957). 

Above  the  condensation  height  hc,  these  equations  were 
augmented  by  one  for  a  wet  adiabat  and  an  explicit  calculation 
of  the  density  defect.  The  numerical  model  was  exercised  using 
the  fire  parameters  in  Table  4  and  the  U.S.  Standard  Atmosphere 
of  1976,  with  a  moisture  burden  given  by  the  midlatitude  average 
profile  (see  Appendix  C,  Fig.  C-2). 

The  results  of  the  computations  are  given  in  Table  4,  where 
the  various  methods  can  be  compared.  The  numerical  method  gives 
results  for  the  wildland  fires  similar  to  those  of  the  empirical 
method,  but  these  do  not  necessarily  represent  reality  either. 
The  (vertical)  rate  of  rise  of  many  of  the  plumes  was  calculated 
as  a  few  cm/s  below  the  point  of  neutral  buoyancy,  so  neglecting 
the  effect  of  the  (horizontal)  10  m/s  windspeed  is  wrong. 

*Projecting  the  burning  area  below  the  surface  to  a  "virtual 
origin"  as  suggested  in  Morton  (1957)  gives  an  origin  farther 
below  the  surface  than  the  plume  height. 


TABLE  4.  PLUME  RISE  CALCULATIONS  FOR  THE  FIRES  OF  TABLE 
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Note  that  the  methods  differ  markedly  in  their  predictions  of 
plume  rise  for  the  Tacoa  oil  fire,  for  which  the  empirical 
method  of  Chandler  et  al.  was  not  designed. 

The  point  here  is  that  the  misapplication  of  a  theory  may 
sometimes  yield  a  believable  result,  but  it  may  not  represent 
reality.  Thus,  when  such  a  theory  is  applied  to  a  situation 
with  which  the  user  is  unfamiliar,  meaningless  results  may  be 
obtained  and  yet  go  unrecognized. 


2.3  A  QUALITATIVE  DISCUSSION  OF  THE  FORMATION,  SPREADING,  AND 
DISAPPEARANCE  OF  A  WATER  CLOUD 


The  discussion  of  Section  2.2  points  out  that  a  large  fire 
creates  a  significant  combustion  plume  which  takes  a  substantial 
amount  of  water  vapor  to  stabilization  at  height  hs .  A  qualita¬ 
tive  discussion  of  the  formation,  spreading,  and  disappearance 
of  a  water  cloud  may  be  based  on  the  following  assumptions. 

(a)  During  plume  rise,  1  to  10  kg  H20/kg  fuel  is  entrained 
and  brought  to  stabilization. 

(b)  Stabilization  occurs  in  the  middle  troposphere,  which 
at  midlatitudes  is  assumed  to  have  50  percent  relative 
humidity  (see,  e.g.,  Appendix  C,  Fig.  C-2). 

(c)  For  a  point  source  of  tracer,  if  one  describes  mean 
wind  dispersion  in  terms  of  turbulence,  the  horizontal 
mean  cloud  width  is  given  in  Fig.  4  as  a  function  of 
travel  time,  where  the  range  of  curves  I-III-IV  is  a 
measure  of  the  variability  of  atmospheric  turbulence. 
We  assume  a  vertical  cloud  width  az  =  0.5  km  and 

for  horizontal  cloud  width  use  the  values  of  Oy 
from  Fig.  4  computed  at  t  =  1  hour  for  curves  I 
(minimum).  III  (mean),  and  IV  (maximum)  which  are, 
respectively,  0.65,  7,  and  60  km.  The  cloud  volume 
for  one  hour  of  burning  is 

V  =  a  2  a  .  (5) 


■r.V.v . 


.  .  »  •  »  *  *  . . «  •  _>»  r 
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HORIZONTAL  MEAN  CLOUD  HALF  WIDTH. 
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FIGURE  4.  Horizontal  dispersion  as  a  function  o£  travel  time. 

The  mean  curves  and  bounds  apply  to  atmospheric 
altitudes  up  to  20  to  25  km.  Earlier  estimate  of 
Hage  (1964),  Hage  et  al.  (1966)  indicated  slower 
spreading  than  do  the  more  recent  data.  Thus,  at 
present,  the  mean  curve  III  with  bounds  I  and  IV 
is  recommended  as  against  Hage's  mean  curve  II  and 
bounds  I  and  III. 


Source:  Bauer  (1983). 
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Now  for  conditions  A-D  above,  Table  5  lists  whether  or  not 
water  clouds  are  formed,  and  if  they  are,  how  long  they  remain 
before  expansion  takes  them  into  a  volume  which  is  no  longer 
saturated  so  that  the  cloud  evaporates.  We  see  that  for  violent 
atmospheric  turbulence  (curve  IV)  the  cloud  expansion  is  so 
rapid  that  condensation  never  occurs,  while  for  minimal  atmo¬ 
spheric  turbulence  (curve  I)  clouds  are  generally  formed,  and 
remain  for  times  of  1  to  10  hours. 

It  is  clear  that  the  assumptions  A,  B,  and  C  are  only  sche¬ 
matic,  but  they  do  lead  to  numbers  that  are  not  unreasonable. 
Referring  back  to  Fig.  2,  we  see  that  the  capping  cumulus  has  a 
maximum  dimension  of  order  20  km  which,  from  Fig.  4  would  sug¬ 
gest  a  mean  lifetime  of  perhaps  two  hours  for  curve  III  or 
"average"  atmospheric  turbulence.  This  is  certainly  not  incon¬ 
sistent  with  the  behavior  of  actual  cumulus  clouds  which  occur 
under  relatively  unstable  conditions,  i.e.,  between  curves  III 
and  IV  of  Fig.  4. 

2.4  LATE-TIME  BEHAVIOR  OF  SMOKE 

The  present  discussion  deals  with  the  generation  of  the 
smoke  cloud,  including  the  formation  and  disappearance  of  a  cap¬ 
ping  cumulus  in  its  first  few  hours.  What  happens  after  that? 

The  (black  or  gray)  smoke  plume  left  behind  after  the  water 
cloud  has  evaporated  is  advected  within  an  "air  parcel".  The 
motion  of  this  parcel  is  approximately  isentropic  but  the  cloud 
absorbs  solar  energy  during  the  daytime  and  in  so  doing  gains 
some  buoyancy,  so  that  it  rises  more  than  it  otherwise  would 
(see,  e.g.,  Bauer,  1984,  Appendix  E).  (Note  that  the  cooling 
due  to  thermal  emission  in  the  IR  will  be  relatively  small,  as 
is  the  greenhouse  effect  due  to  the  smoke  (see,  e.g.,  Bauer, 
1984,  Appendix  C,  McCartor  et  al.,  1985).  Thus,  the  air  parcel 
will  gain  perhaps  90  percent  of  the  absorbed  solar  energy,  lead¬ 
ing  to  a  net  increase  in  its  potential  temperature.) 
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Eventually  the  smoke  must  be  scavenged  from  the  atmosphere, 
and  there  are  at  least  three  distinct  mechanisms  with  possibly 
different  time  scales: 

1.  Immediate  rainout  leading  to  "black  rain".  This  has 
been  observed  (e.g.,  in  Dresden  and  Hiroshima  in  World 
War  II  bombing  raids,  and  also  in  some  forest  fires). 
Calculations  (e.g.,  Molenkamp,  1983)  agree  with  obser¬ 
vations  which  suggest  that  this  process  only  occurs 
some  10  percent  of  the  time  and  so  removes  less  than 
this  fraction  of  smoke  on  the  average.  Certainly  the 
present  discussion  does  not  offer  a  mechanism  for 
immediate  rainout,  though  one  can  imagine  that  it  might 
occur  if  a  fairly  fresh  smoke  plume  encounters  a  parcel 
of  moist,  unstable  air. 

2.  Wet  or  dry  transport  to  the  surface  and  subsequent 
removal.  The  time  scale  for  this  is  one  to  three 
weeks  for  most  non-reactive  tracers  in  the  troposphere 
(see,  e.g.,  Pruppacher  and  Klett,  1978),  somewhat 

less  for  reactive  materials  such  as  sulfuric  and  nitric 
acids.  The  aged  smoke  particles  will  presumably 
disappear  as  fast  as  the  average  "non-reactive  tracer" 
since  they  will  not  be  as  hydrophobic  as  when  they  are 
fresh . 

3.  Heterogeneous  chemical  removal,  e.g.,  by  reaction  with 
ozone  or  other  reactive  atmospheric  material.  Prelim¬ 
inary  observations  (deHaas  et  al.,  1986,  and  Golden, 
1986)  suggest  characteristic  times  of  one  to  three 
months  (with  uncertainty  by  perhaps  an  order  of  magni¬ 
tude  up  or  down)  for  soot  removal. 

2.5  OPTICAL  PROPERTIES  OF  SMOKE 

As  has  been  pointed  out  above,  the  mass  of  a  capping  cum¬ 
ulus  or  other  cloud  of  entrained  and  condensed  water  is  much 
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greater  than  that  of  a  late-time  black  or  white/gray  smoke 
plume,  because  the  mass  of  water  from  both  combustion  and 
entrainment  is  several  orders  of  magnitude  greater  than  the 
mass  of  smoke.  Note,  however,  that  this  condensed  water  will 
evaporate  rather  rapidly  (see,  e.g.,  Section  2.4)  so  that  after 
some  hours  it  will  all  be  gone. 

Thus,  the  climatic  impact  of  the  "Nuclear  Winter"  smoke 
cloud  is  largely  that  of  the  smoke,  not  of  the  condensed  water. 
(Indeed,  it  is  now  considered  that  the  net  climatic  impact  of 
water  clouds  is  rather  small,  in  that  the  cooling  due  to  the 
scattering  of  sunlight  by  the  cloud  roughly  balances  the  green¬ 
house  warming  due  to  the  absorption  of  earthshine  (see,  e.g.. 
Cess  et  al . ,  1982 ) . 

Here  we  shall  review  the  optical  properties  of  water  and 
smoke  particles  in  both  visible  and  infrared,  distinguishing 
between  smoke  from  oil  and  wood  fires.  For  reference.  Table  6 
contains  the  complex  refractive  index  of  water  and  soot,  as 
well  as  of  some  other  atmospheric  aerosols.  Note  that  in  the 
visible,  soot  is  much  more  absorptive  than  water,  whereas  in 
the  IR,  especially  for  X  Z  10  um,  most  materials  have  similar 
absorptivity  N2. 

The  smoke  from  oil  fires  looks  black  and  is  here  assumed 
to  be  made  up  of  soot,  which  is  typically  considered  to  be 
made  up  of  very  small  particles  (>  0.01  um  minimum  dimension) 
which  tend  to  clump  together.  Here  these  soot  particles  are 
treated  as  spheres  of  radius  a,  which  are  characterized  by  the 
Mie  parameter  q: 

q  =  2ir  a/X  (  6  ) 

for  light  of  wavelength  X.  Here  q  <<  1,  i.e.,  we  are  in  the 
Rayleigh  region  in  which  there  is  very  little  scattering,  while 
the  absorption  cross  section  is  proportional  to  the  volume  of 
the  particle  divided  by  the  wavelength  of  light.  Explicitly 
(cf.  Van  de  Hulst,  1957), 
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TABLE  6.  REFRACTIVE  INDEX  FOR  VARIOUS 
AEROSOL  COMPONENTS 
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=  ( 32/27 )  n  a2  q4  (N1  -  1  ) 


=  ( 8/3 )  *  a2  qN  , 


where  N^  and  N2  are,  respectively,  the  real  and  imaginary  parts 
of  the  complex  refractive  index. 

The  smoke  from  oil  fires  is  quite  variable  in  both  amount 
and  characteristics  depending  on  the  characteristics  of  the 
fires  (cf,  (e.g.,  presentations  at  the  20th  Symposium  (  Inter- 
nat lona 1 )  on  Combustion,  Ann  Arbor,  MI,  August  1984).  It 
often  contains  a  significant  number  of  large  particles,  which 
scatter  as  much  as  they  absorb.  Note  however: 

(a)  the  assumption  that  much  of  the  absorption  comes 
from  small  particles  for  which  Eq.  (7)  holds  may 
still  be  useful  if  the  electric  field  vector  of 
the  incident  radiation  is  parallel  to  the  smallest 
dimension  of  the  soot  particles. 

(b)  while  a  capping  cloud  exists,  one  will  tend  to  see 
this  rather  than  the  smoke  cloud  because  of  its  much 
larger  mass  (see,  e.g..  Fig.  2). 

Fraser  et  al.  (1986)  have  analyzed  visible  satellite 
imagery  (from  GOES)  of  the  July  1982  Yukon  forest  fires  over 
the  Great  Lakes,  over  Chesapeake  Bay,  and  over  the  North  Atlan¬ 
tic.  Because  of  the  long  transit  time  ( 1  to  3  days)  this  smoke 
has  been  aged  so  that  any  capping  cumulus  clouds  will  presum¬ 
ably  have  evaporated.  Fraser  et  al.  (1986)  find: 

(a)  maximum  optical  thickness  ~  1.6  -  2  in  vertical 
viewing 

(b)  single  scattering  albedo  w ,  defined  as 


w  =  °sca  /  (°sca  +  aabs' 


in  the  range  0.92  -  0.98,  i.e.,  extinction  is  due 

mainly  to  scattering  rather  than  to  absorption 


(c)  mean  particle  size  a  ~  0.3  um,  so  that  the  Mie 
parameter  q  =  2ira/X~2  in  the  visible,  which 
indicates  that  scattering  can  be  significant  (if 
q  <<  1,  ogca  a  q4,  which  is  very  small). 

It  must  be  stressed  that  the  results  of  Fraser  et  al .  have 
been  obtained  by  an  as-yet-unvalidated  method.  Nevertheless, 
they  do  suggest  that  a  forest  fire  of  the  Yukon  (boreal)  type 
produces  smoke  which  is  largely  scattering  rather  than  absorp¬ 
tive  in  the  visible — and  GOES  imagery  of  the  visible  plume  over 
the  North  Atlantic  on  2  August  1982  is  consistent  with  that 
result . 

One  question  which  must  be  addressed  is  the  following: 
can  one  tell  "by  eye"  whether  a  given  cloud  is  largely  absorp- 
tive/( "black" )  or  scattering/( "white" ) ?  The  answer  is  somewhat 
subtle.  Table  7  gives  the  effective  hemispheric  reflectivity 
for  optically  thick  clouds  as  a  function  of  single  scattering 
albedo  or  the  ratio  of  absorption  to  scattering  cross  section. 
The  reflectivity  clearly  depends  on  the  angle  of  illumination: 
for  near-normal  illumination  the  reflectance  is  smaller  than 
for  near-tangential  (glancing  angle)  illumination.  However, 
even  for  crat>s/asca  ~  0.15  to  0.20  the  hemispherical  reflectance 
for  near-normal  illumination  is  of  order  0.1  (for  tangential 
illumination  it  would  be  ~  0.25).  Thus,  even  if  a  large  part 
of  the  extinction  is  due  to  scattering,  the  plume  or  cloud  may 
still  look  dark  gray. 

This  raises  a  major  unknown — what  is  smoke?  Is  it  largely 
an  absorbing  medium  made  up  of  graphitic  soot — particles  so 
small  that  they  do  not  scatter  much,  i.e.,  q  =  2:ra/X  <  1,  e.g., 
a  ~  0.01  um — or  in  part  larger  scattering  particles  for  which 
q  £  1,  i.e.,  a  2  0.09  um  for  X  <,  0.55  um.  The  question  is 
significant  because  if  q  <<  1,  asca  ~  q4  n  a^  -*•  0  so  that 
aext  =  asca  +  aabs  ~  aabs  £rom  E(3*  (7),  and  the  effective 
optical  thickness  Xe  of  a  cloud  of  geometrical  thickness  L  and 
mass  loading  M  (g/m^)  is 
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RIC  REFLECTANCE  FOR  OPTICALLY  THICK  CLOUDS  AS  A  FUNCTION  OF 


(  3 tt / 2  )  N2  ml  /  Xp 


for 


q  <<  1 


(10) 


where  n  =  number  density  of  particles  of  density  =  M/ ( 4 n p a -V 3 ) . 
By  contrast,  if  q  ^  1,  °ext  ~  C1T“  (c  ~  2)  so  that 

Xe  =  3cML/4ap  for  q  ^  1  (11) 

so  that  if  q  2  1,  larger  particles  produce  a  significantly 
smaller  optical  thickness  and  thus  extinction,  cooling,  etc., 
than  do  smaller  particles. 

Evidently,  either  absorbing  or  scattering  particles — or 
any  combination  thereof — can  produce  surface  cooling,  but  the 
numerical  magnitude  of  the  effect  may  depend  critically  on  some 
details  of  what  the  smoke  actually  is. 

2.6  CLIMATIC  EFFECTS  OF  SMOKE 

(a)  The  climatic  impact  of  water/ice  clouds  has  been 
studied  by  many  people,  in  particular  by  Cess  et  al. 
(1982).  The  key  question  is  the  relative  importance 
of  cooling  due  to  reduction  in  incoming  solar  radia¬ 
tion  ("albedo  effect")  and  the  reduction  in  outgoing 
earthshine  ("IR  radiation").  The  conclusion  is  that 
these  two  relatively  large  effects  almost  balance,  so 
that  there  is  either  a  very  small  net  heating  or  a 
very  small  net  cooling — but  different  data  give  dif¬ 
ferent  signs  of  the  net  effect.  Whether  the  clouds 
are  high  or  low  can  also  be  important  in  this  regard. 

(b)  As  far  as  the  absorption  from  soot  is  concerned  we 
assume  that  Mie  parameter  q  of  Eq.  (6)  is  very  much 
less  than  one  in  the  visible  as  well  as  in  the  LWIR. 
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Thus,  the  optical  thickness  in  the  LWIR,  X(LW), 
is  very  much  smaller  than  that  in  the  visible,  X(VIS), 
because 


X  ( LW ) 
X( VIS) 


q(LW) 
q(  VIS 


N2(LW) 

X  N2( VIS) 


where  N2  =  imaginary  part  of  complex  refractive  index; 
typically  N2 ( LW) /N2 ( VIS )  ~  3  to  10  whereas  q(LW)/q(VIS) 
=  0.05.  Accordingly,  the  primary  effect  of  soot/ 
smoke  particles  is  cooling  due  to  extinction  of  solar 
radiation  rather  than  greenhouse  warming. 

The  climatic  effect  of  grey/white  smoke  is  harder  to 
estimate.  In  the  visible,  solar  radiation  will  be 
extincted  by  scattering  rather  than  by  absorption, 
while  the  greenhouse  warming  in  the  LWIR  will  probably 
be  small.  However,  the  quantitative  estimate  of  net 
cooling  due  to  gray/whit.e  clouds  will  presumably 
depend  on  details  of  particle  size  and  absorptivity 
(see,  e.g.,  Bartky  and  Bauer,  1966,  and  Bartky,  1968). 


3.  DISCUSSION 


This  analysis  demonstrates  how  the  initial  plume  rise  to 
stabilization  depends  on  the  intensity  of  the  fire  and  also  on 
atmospheric  conditions:  a  moist  lower  atmosphere  leads  to  an 
enhanced  plume  rise  and  sometimes  to  the  formation  of  a  capping 
cumulus  cloud  which  is,  however,  relatively  short-lived  (to 
several  hours — see  Fig.  2  or  Table  5) .  Because  of  this  short 
lifetime,  the  climatic  impact  of  the  capping  cumulus  cloud  is 
small . 

However,  after  the  capping  cumulus  or  other  water  cloud 
evaporates  due  to  expansion  into  the  dry  atmosphere,  the  resid¬ 
ual  smoke  plume  continues  to  rise  under  its  solar-induced 
buoyancy  (cf,  e.g.,  Bauer,  1984,  p.  E-4,  especially  Eq.  (4), 
also  Knox  et  al.,  1986,  p.  20),  so  that  the  final  plume  rise 
height  will  tend  to  be  specified  by  its  late-time  behavior 
rather  than  by  its  initial  injection  height.  Note  that  the 
"Nuclear  Winter"  cooling  from  a  high-altitude  smoke  plume  is 
significantly  greater  than  that  from  a  low-altitude  smoke  plume 
(e.g.,  compare  the  results  of  Ramaswamy  and  Kiehl,  1985,  Fig. 
17),  which  after  20  days  show  a  cooling  of  33°C  for  the  high- 
altitude  (constant  density)  profile  D  as  compared  to  23°C  cool¬ 
ing  for  the  lower-altitude  (constant  mixing  ratio)  profile  M. 

In  addition,  the  residence  time  of  material  in  the  atmosphere 
is  greater  for  higher  altitudes,  so  that — other  things  being 
equal  —  the  higher-altitude  smoke  distribution  D  will  remain 
longer  than  the  lower-altitude  profile  M. 

By  looking  at  Fig.  2  or  at  any  instantaneous  photograph  of 
a  smoke  plume,  one  notes  that  the  smoke  on  a  given  occasion  is 


. .  w  ,  ■  'V' 


not  uniformly  distributed  in  a  horizontal  plane,  although  the 
distribution  averaged  over  many  different  smoke  injections  will 
be  uniform.  It  can  be  shown  that  if  a  given  amount  of  smoke  is 
distributed  uniformly,  its  total  extinction  of  the  underlying 
surface  will  have  a  maximum  value.  If  the  distribution  is  non- 
uniform  ("lumpy"),  more  radiant  energy  will  penetate  the  cloud 
to  give  less  net  cooling.  Appendix  D  gives  a  simple  example 
for  a  plume  of  optical  thickness  X  and  transmission  e“x  (i.e., 
absorption  only):  if  X  ~  3  to  10,  the  effective  extinction  can 
be  reduced  by  a  factor  of  1.5  to  8  if  the  thick  and  thin  plumes 
differ  in  optical  thickness  by  a  factor  of  2  to  16.  This  is 
clearly  a  very  simple  example:  a  more  refined  calculation  will 
give  different  numerical  results,  but  the  most  effective  way 
to  solve  this  problem  is  experimentally,  by  studying  the  trans¬ 
mission  through  actual  smoke  plumes. 

An  additional  factor  of  importance  is  provided  by  the  opti¬ 
cal  properties  of  the  smoke,  in  particular  by  the  ratio  of 
absorption  to  scattering,  where  absorbing  (black,  oil)  smoke 
tends  to  give  less  transmission  and  thus  more  cooling  than  does 
scattering  (white,  wood)  smoke. 

From  the  original  photo  on  which  Fig.  2  is  based,  one  can 
see  that  during  land  clearing  in  Brazil  in  August  1984  the  smoke 
plume  remaining  after  the  capping  cloud  over  the  large  fire 
evaporates  is  white  or  light  gray,  not  black,  i.e.,  largely 
scattering  rather  than  absorbing.  This  is  indeed  characteristic 
of  smoke  from  wood  fires  (forest  or  agricultural).  In  terms  of 
the  single  scattering  albedo 

u  =  CT sea  /  0 ext  •  (12) 


where 


a 


ext 


asca  +  CTabs  • 


(13) 


Fraser  et  al.  (1986)  find  u  =  0.92  -  0.98  for  smoke  from  the 
Yukon  fires  of  July  1982  after  1  to  3  days'  travel,  using  an 
as-yet-unval idatd  method  of  analyzing  visible  GOES  satellite 
data.  Ginsburg  and  Golitsyn  (1986)  quote  w  =  0.9  for  represen¬ 
tative  fires  in  the  USSR. 

By  contrast,  Roessler  and  Faxvog  (1981)  quote  to  ~  0.15  for 
soot,  i.e.,  the  highly  absorptive  smoke  from  oil  fires.  (Such 
low  values  are  supported  qualitatively  from  observations  on  oil 
pool  fires--B.  Zak,  Sandia  National  Laboratory,  private  communi¬ 
cation.  ) 

From  an  analysis  of  the  available  fuel  in  a  "representative" 
city,  Knox  et  al.  (1986,  p.  19)  suggest  that  15  to  20  percent  of 
fuel  consists  of  oil  or  oil-derived  material,  with  the  remainder 
being  wood  or  wood-derived.  If  the  wood-derived  smoke  has  u  ~ 

0.9  to  0.95  and  if  the  oil  has  no  scattering,  i.e.,  w  =  0.0, 
and  all  materials  have  the  same  smoke  emission  index,  this  gives 
w  ~  0.7  to  0.8  for  a  major  urban  fire. 

The  effect  of  this  contribution  of  scattering  to  extinction 
is  to  increase  the  transmission  through  a  cloud  of  given  extinc¬ 
tion  optical  thickness  Xe : 

xe  =  aext  k  ,  (14) 

where  n  =  number  of  particles  per  unit  volume  and  L  =  geometrical 
cloud  thickness.  Table  8  gives  the  transmission  coefficient 
t(X;u)  for  a  plane  wave  incident  normally  on  a  uniform  medium 
as  a  function  of  optical  thickness  X  and  single-scattering  albedo 
The  results  show  that  for  a  given  optical  thickness  X  absorbers 
( to  small)  transmit  less  radiant  energy  than  do  scatterers  (w 
close  to  1),  so  that  the  effective  "Nuclear  Winter"  cooling  due 
to  absorbing  (black,  oil)  smoke  is  greater  than  that  due  to  scat¬ 
tering  (white,  wood)  smoke.  The  calculations  are  given  for  two 
simple  models  A  (isotropic)  and  B  (forward  +  backward  scattering 
only)  to  demonstrate  these  qualitative  characteristics:  the  numer 
ical  values,  of  course,  depend  on  the  particular  model  adopted. 


TABLE  8.  TRANSMITTANCE  IN  VISIBLE  AS  A  FUNCTION  OF  SINGLE  SCATTERING  ALBEDO 
EXTINCTION  OPTICAL  THICKNESS  ( ) ,  FOR  TWO  DIFFERENT  MODELS  A  AND  ! 
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APPLICATION  TO  THE  GLOBAL  EFFECTS/NUCLEAR  WINTER  PROBLEM 

The  preceding  discussion  pointed  out  the  following: 

(a)  For  the  climatic  cooling,  the  initial  plume  rise  may 
well  be  less  important  than  the  long-term  motion  due 
to  solar-induced  buoyancy  and  other  atmospheric 
motions.  (Ramaswamy  and  Kiehl,  1985,  indicate  a 
possible  variability  by  a  factor  of  1.5  between  a 
"constant  density"  and  a  "constant  mixing  ratio" 
model . ) 

(b)  On  any  given  occasion,  the  smoke  plume  will  be  quite 
nonuniform.  The  result  will  be  to  reduce  the  effec¬ 
tive  cooling  by  a  significant  factor  (the  simple 
model  of  Appendix  D  suggests  a  factor  of  1.5  to  8, 
but  this  is  highly  uncertain). 

(c)  Since — for  a  given  optical  depth — absorption  is  much 
more  effective  than  scattering  in  reducing  transmis¬ 
sion  and  thus  producing  cooling,  the  late-time  optical 
properties  of  the  smoke  are  crucial.  (Table  8  suggests 
that  a  variability  by  a  factor  of  2  to  5  is  likely  for 
transmission  between  different  models  that  have  been 
used .  ) 

How  can  these  effects  be  studied? 

(a)  Perhaps  the  best  example  of  long-term,  solar-induced 
buoyancy  is  provided  by  Arctic  Haze.  As  an  example, 
the  AGASP  program  (see,  e.g.,  Schnell,  1984)  demon¬ 
strates  that  the  plume  of  the  Norilsk  copper-nickel 
smelter  could,  on  occasion,  be  followed  for  10  days 
from  the  surface  source  in  Siberia  to  layers  near  3 


1984;  Ottar  and  Pacyna,  1984;  and  Radke  et  al.,  1984) 

(b)  Smoke  nonuniformity  can  be  studied  by  quantitative 
photographic  or  video  techniques  by  looking  at  plumes 
the  plumes  from  larger  sources  can  normally  be 
studied  for  longer  times  than  those  from  smaller 
sources,  depending  on  atmospheric  conditions. 

(c)  The  effective  climatic  cooling  will  be  determined  by 
late-time  optical  properties  of  smoke,  and  this 
requires  large  sources  which  can  be  followed  for  long 
times.  For  wood/forest  fires,  the  best  (if  unvalida¬ 
ted)  method  is  that  used  by  Fraser  et  al.  (1986)  for 
from  1  to  3  days.  Sooty  smoke  from  oil  fires  may 
react  with  ozone  at  very  long  times  (months — cf.  e.g. 
de  Haas  et  al.,  1986,  and  Golden,  1986);  the  best 
approach  right  now  is  to  use  the  plume  from  large  oil 
pool  fires  at  Sandia  National  Laboratory,  which  can 
be  tracked  with  an  aircraft  for  several  hours  (Zak  et 
al.,  1986). 


5. 


RECOMMENDATIONS 


1.  Use  Arctic  Haze  data  to  study  long-term  plume  rise  by  solar 
induced  buoyancy  and  transport. 

2.  Use  data  from  large  fires  to  study  the  evolution  of  the 
late-time  plume  nonuniformity  which  will  affect  the  trans¬ 
mittance  of  solar  radiation  and  thus  the  long-term  climatic 
cool i ng . 

3.  Validate  the  method  of  Fraser  et  al.  (1986)  for  studying 
the  late-time  (1-  to  3-day)  optical  properties  of  smoke 
from  wood  fires  using  weather  satellite  imagery. 

4.  Use  the  large  Sandia  National  Laboratory  oil  pool  fires  as 
a  source  to  study  late-time  (1-  to  6-hour)  optical  proper¬ 
ties  of  sooty  oil  smoke  in  a  dry  atmosphere. 
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TASK  ORDER 

NO.  MDA903  84  C  0031:  T-4-237 

TITLE:  Atmospheric  Injection  and  Scavenging  of 
Smoke  Particles  from  Fires 


1.  This  task  is  for  work  to  be  performed  by  the  Institute  for 
Defense  Analyses  for  the  Defense  Nuclear  Agency. 

2.  BACKGROUND: 

The  current  concern  of  major  climatic  cooling  arising  from  a 
nuclear  exchange  is  based  largely  on  the  scenario  that  nuclear 
detonations  will  ignite  large  fires.  These  fires  produce  a  large 
quantity  of  fine  smoke  particles,  which  -  when  carried 
sufficiently  high  in  the  atmosphere  -  are  claimed  to  produce 
cooling  at  the  earth's  surface.  The  altitude  of  the  smoke  cloud 
is  alleged  to  determine  whether  cooling  or  warming  occurs,  while 
the  rate  of  removal  of  the  smoke  particles  will  determine  the 
duration  of  the  climatic  effects.  An  evaluation  of  the  quantitie 
of  smoke  particles,  the  rise  height,  and  the  scavenging  rate  is 
needed  as  a  first  step  in  determining  the  credibility  of  this 
scenario.  To  do  so,  data  on  different  types  of  fires  (urban, 
forest,  oil,  etc.)  and  on  removal  mechanisms  and  rates  will  be 
examined. 

3.  OBJECTIVE  AND  SCOPE: 

To  evaluate  the  magnitude  of  particulate  injection,  the  rise 
height,  and  the  rate  of  removal  of  smoke  particles  put  into  the 
atmosphere  as  a  consequence  of  large-scale  fires  as  posited  to 
result  from  a  massive  nuclear  exchange.  The  work  will  review 
current  knowledge  and  uncertainties  and  make  recommendations  on 
how  the  uncertainties  may  be  reduced. 

4.  SCHEDULE: 

This  effort  will  begin  1  January  1984.  A  draft  document  will 
be  delivered  by  30  September  1984,  and  a  final  version  90  days 
thereafter.  Memoranda  and  briefings  will  be  given  as  significant 
portions  of  the  effort  are  completed. 


ACSCAKCH  ANO 
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.  TECHNICAL  COGNIZANCE: 

Technical  cognizance  for  this  task  is  assigned  to  the 
Director ,  Radiation  Directorate,  Defense  Nuclear  Agency, 
Attn:  RAAE. 

6.  FUNDING: 

The  expenditure  of  $74,000  of  FY  1984  funds  is 
authorized  for  this  task. 

7.  SPECIFIC  ADMINISTRATIVE  INSTRUCTIONS: 

a.  If  at  any  time  during  the  course  of  this  task  IDA 
identifies  the  need  for  changes  in  this  task,  such  as 
additional  resources,  schedule  modification,  changes  to 
emphasis  of  effort  or  scope,  etc.,  as  set  forth  in  the  above 
paragraphs,  a  report  with  appropriate  recommendations  will 
be  submitted  in  accordance  with  the  terms  of  the  IDA/WSEG 
Memorandum  of  Understanding  of  12  March  1975  (and  its 
successor)  as  applicable  to  the  Director,  DoD-IDA  Management 
Office,  OUSDRE,  with  a  copy  to  the  sponsor  or  his  project 
officer,  as  appropriate.  Changes  in  this  task  will  be  made 
only  with  the  approval  of  appropriate  cognizant  DoD 
officials. 

b.  This  task  will  be  conducted  under  Industrial  Security 
Procedures  in  the  IDA  area.  If  certain  portions  of  the  task 
require  the  use  of  sensitive  information  which  must  be 
controlled  under  military  security,  the  DoD-IDA  Management 
Office  will  provide  supervised  working  areas  in  which  work, 
will  be  performed  under  military  security  control. 

c.  The  termination  date,  the  date  after  which  no  further 
costs  will  be  incurred  against  the  contract,  for  this  task 
is  7  January  1985;  unless  changed  by  a  written  amendment  to 
this  task  order. 

d.  A  "need  to  know"  is  hereby  established  in  connection 
with  this  task.  Specifically,  authority  is  granted  to  IDA 
to  request,  in  my  name,  classified  and  unclassified 
documents  and  publications  needed  to  accomplish  this  task 
and  to  verify,  using  my  name,  the  "need  to  know"  with 
respect  to  access  to  classified  information  needed  to 
complete  this  task.  Further,  IDA  authority  to  request 


OFFICE  OF  the  under  secretary  of  DEFENSE 
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12  April  1985 

HCSCAftCH  AND 
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TASK  ORDER 

NO.  MDA903  84  C  0031:  T-4-237 

AMENDMENT  NO.  1 

TITLE:  Atmospheric  Injection  and  Scavenging  of  Smoke 
Particles  from  Fires 


1.  This  task  amendment  is  to  extend  work  initiated  under  this 
task.  Paragraphs  2,  3,  4,  6,  and  8  are  changed  as  follows: 

2.  BACKGROUND:  Add: 

"Task  T-4-237  reviewed  how  much  smoke  is  generated  per 
unit  area  by  urban  and  forest  fires,  and  some  critical  aspects  of 
atmospheric  injection.  It  indicated  the  significance  of 
meteorological  variability,  and  how  a  study  of  forest  fires  can 
help  to  indicate  what  happens,  in  particular  the  effective  smoke 
injection  altitude  and  rate  of  disappearance. 

Overall,  the  "Nuclear  Winter"  scenario  has  not  been 
shown  to  be  physically  impossible,  and  thus  it  is  appropriate  to 
continue  investigating  all  relevant  aspects  of  the  scenario  with 
emphasis  on  actual  physical  phenomena  which  simulate  different 
aspects  of  Nuclear  Winter.  One  phase  which  has  not  been 
investigated  in  detail  is  how  smoke  is  dispersed  from  its 
originating  fires,  near  the  ground,  to  the  globally  uniform 
distribution  throughout  the  troposphere  which  is  used  in  climate 
models  as  a  basis  for  the  predictions  of  a  very  large  cooling. 

It  is  known  that  large-scale  smoke  and  smog  plumes  travel  over 
large  distances,  both  horizontally  and  vertically,  and  here  a 
study  of  forest  fires  and  of  arctic  haze  provides  an  opportunity 
to  study  the  long-range  transport  and  disappearance  of  the 
smoke,  in  particular  soot. 

3.  OBJECTIVE  AND  SCOPE:  Add: 

"Work  under  Amendment  1  will  use  available  data  on  forest 
fires,  arctic  haze,  and  other  relevant  phenomena  which  provide 
partial  simulations  for  the  Nuclear  Winter  scenario  with  the 
objectives  of: 

a.  Bounding  the  effects 

b.  Providing  a  baseline  for  computer  models  of  cloud 
spreading 

c.  Providing  input  information  for  climate  models 


Ongoing  experimental  and  analytical  work  at  NASA,  NOA. 
etc.,  will  be  used  and  extended  as  feasible,  and  recoramendati< 
on  further  work  will  be  made.” 

4.  SCHEDULE:  Add: 

"This  effort  under  Amendment  1  will  begin  on  1  Decembei 
1984.  A  draft  document  will  be  delivered  by  30  September  198! 
and  a  draft  final  version  90  days  thereafter.  Memoranda  and 
briefings  will  be  given  as  significant  portions  of  the  effort 
are  completed.” 

6.  FUNDING:  Add: 

"The  expenditure  of  $125,000  of  FY  1985  funds  is 
authorized  to  continue  work  on  this  task,  as  amended.” 

8.  SPECIFIC  ADMINISTRATIVE  INSTRUCTIONS:  Replace 
subparagraph  c.  with  the  following:  ” 

"c.  The  termination  date,  the  date  after  which  no  furtt 
costs  will  be  incurred  against  the  contract,  for  this  task  is 
December  1985.” 


ACCEPTED : 


T.  L.  RICKETTS 
Colonel  USA 
Director 

DoD-IDA  Management 


A.  J.  GOOD  POSTER 
General,  U.  S.  Army  (Ret.) 

President,  Institute  for  Defense  Analyses 


Office 
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May  13.  1985 


INSTITUTE  FOR  DEFENSE  ANALYSIS 

1801  S  Beauregard  Si  reel .  Alexandria.  Virginia  22*11  •  Telephone  (70.*i  84^  2(*n. 
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January  14,  1986 


MEMORANDUM  FOR  THE  RECORD 

SUBJECT:  Task  Order  T-4-237,  Amendment  1,  Atmospheric 

Injection  and  Smoke  Particles  from  Fires 


In  accordance  with  previous  discussions,  it  is  agreed  that 
the  final  work  on  this  Amendment,  on  black  and  white  clouds  from 
oil  fires  and  implications  for  nuclear  winter,  will  be  delayed 
at  no  cost  to  the  government  until  1  August  1986.  At  that  time, 
a  draft  report  will  be  furnished  with  a  proposed  distribution 
list.  DNA  ( RAAE )  will  return  comments  on  the  report,  including 
approval  of  classification  and  approved  distribution  list, 
within  three  months  thereafter.  The  report  will  be  published 
by  IDA  in  final  form  three  months  after  receipt  of  sponsor's 
comments.  Printing  and  distribution  will  be  made  when  an 
approved  distribution  list  is  furnished  to  IDA,  using  resources 
provided  by  the  task. 


Robert  E.  Roberts 
Director 

Science  and  Technology  Division 
Institute  for  Defense  Analyses 


Leon  A.  Wittwer 
Program  Manager 
Atmospheric  Effects  Division 
Defense  Nuclear  Agency 


APPENDIX  B 


SUMMARY  OF  IDA  MEMORANDUM  REPORT  M-116 

This  paper  suggests  a  reason  and  a  means  to  study  large- 
scale  natural  and  urban/industrial  fires.  Very  large  fires — 
both  forest/agricultural  and  urban/industrial--provide  a 
significant,  if  partial  simulation  for  the  urban  fire  scenario 
depicted  as  playing  a  role  in  the  Nuclear  Winter  hypothesis. 

Only  the  largest  naturally  occurring  fires  provide  the  large- 
scale  atmospheric  smoke  injections  and  perturbations  that  can 
be  traced  over  long  times  (hours  to  days).  Such  fires  occur 
only  very  rarely  (fortunately),  and  thus  they  have  to  be  treated 
as  experiments  of  opportunity,  i.e.,  to  be  studied  when  they 
occur,  using  observations  from  existing  operational  sensors. 
Large-scale/long-time  current  phenomena  are  best  studied  using 
earth/atmosphere  observation  satellites.  Weather  satellites 
provide  global  coverage  with  frequently  updated  imagery,  but  of 
relatively  poor  resolution  (>  1  km)*. 

To  extract  useful  information  from  such  experiments  of 
opportunity  calls  for  pre-planning.  As  an  initial  phase,  a 
workshop  involving  experts  within  DOD,  NOAA  and  NASA  was  held 
at  IDA  on  1  May  1985.  Its  product  -  a  concept  document  on 
experiments  of  opportunity  for  Nuclear  Winter  -  is  reproduced 
here  as  Part  I.  It  concluded  that  when  an  appropriate  fire  is 
reported,  an  ad  hoc  panel  should  be  alerted  to  arrange  in  real 
time  for  all  relevant  satellite  and  other  data  to  be  taken  and 
retained  for  later  analysis.  The  panel  and  the  procedure  for 


*"Regular"  earth  observation  satellites  such  as  LANDSAT  and 
SPOT  provide  much  higher  resolution  (20  to  80  in),  but  the 
typical  18-day  revisit  time  for  any  point  on  the  globe  means 
that  there  is  unlikely  to  be  imagery  of  the  (rare)  events 
under  consideration  here. 


alerting  the  panel  and  arranging  for  data  collection  must  be 
established  in  advance,  and  checked  out  with  a  test  exercise. 

In  the  Nuclear  Winter  scenario,  most  smoke  is  now  consid¬ 
ered  to  be  due  to  urban/industrial  fires  rather  than  to  forest 


fires  because  of  the  total  mass  of  material  burned  and  thus 


the  smoke  generated.  Cities  and  industrial  complexes  have  high 
value  and  are  designed  to  minimize  the  risk  of  devastating 
fires.  By  comparison,  forests  have  a  relatively  low  value  per 
unit  area.  Thus,  it  is  not  surprising  that  large  forest  and 
agricultural  fires  are  larger  than  even  the  largest  industrial 
fires.  Table  B-l  shows  some  representative  large  fires  that 
occurred  in  the  1970-1984  time  frame.  The  industrial/urban 
fires  are  chosen  from  among  the  20-odd  largest  such  in  the 
Americas  over  this  time  period;  slash-and-burn  occurs  each  year 
just  before  the  planting  season  in  Latin  America,  while  a  10^ 
hectare  (2.5  million  acre)  boreal  forest  fire  occurs  in  Yukon, 
Alaska  or  Siberia  at  least  every  10  to  20  years. 

From  Table  B-l  the  question  arises,  how  small  a  fire  can 
be  observed  using  weather  satellite  imagery.  Part  II  of  IDA 
Memorandum  Report  M-116  entitled  "Remote  Sensing  of  Smoke 
Plumes  in  the  Visible:  How  Small  a  Plume/Fire  Can  One  Expect 
to  Detect?"  addresses  this  question.  After  a  parametric  dis¬ 
cussion  of  how  a  fire  is  detected  using  curren.  satellite 
imagery,  the  example  of  the  large  oil  fire  that  occurred  on 
20-22  December  1982  at  Tacoa ,  Venezuela  is  discussed.  Figure 
B-l  shows  the  early  morning  (7.08  a.m.  LT)  high  resolution  (1 
km)  NOAA-6  image  on  20  December  1982;  the  plume  over  Caracas, 
Venezuela  is  due  to  this  fire.  Figure  B-2  shows  the  afternoon 
(2.42  p.m.)  low  resolution  (2  km)  NOAA-7  image  on  the  same  day. 
Here  the  smoke  plume  cannot  be  seen,  possibly  because  it  has 
dissipated  during  the  seven-hour  interval,  but  mainly  because 
of  the  typically  heavy  cumulus  cloud  buildup  in  the  afternoon, 
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N3AA-6  LS  Rev.  18092,  20  Dec  1982,  1208GMT 


Visual  Band,  1.0  km  resolution. 


(Received  from  CIRES/NOAA  National  Snow  & 


Ice  Data  Center/CMSP — presumably  from  Panama 


City  Ground  Station.) 


The  smoke  plume  over  Caracas,  Venezuela, 


presumably  comes  frcm  the  Tacoa  oil  fire. 


Its  dimensions  are  120  km  N/S  (along-wind) 
and  70  km  E/W  (cross-wind),  corresponding 
to  several  hours  of  dispersion. 


m 


FIGURE  B-l. 


Visible  Satellite  Imagery  from  NOAA-6,  7:08  a.m. 
LT  on  20  December  1982  (1  km  resolution)  showing 
the  smoke  plume  from  the  Tacoa,  Venezuela  oil 
fire  over  Caracas. 
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supplemented  by  the  effects  of  lower  resolution  and  a  less 
favorable  sun-cloud-satellite  scattering  angle. 

One  may  draw  the  following  conclusions  from  this  study: 

1.  Large  forest  fires  can  provide  a  useful,  if  partial, 

simulation  of  some  aspects  of  Nuclear  Winter,  in 
particular:  smoke  emission,  plume  spreading,  non¬ 

uniformity  and  dispersion/disappearance. 

2.  Weather  satellites  provide  a  useful  and  inexpensive 
if  incomplete  way  of  studying  the  time/space  develop¬ 
ment  of  a  series  of  very  large  smoke  plumes  from 
forest/agricultural  or  industrial/urban  fires. 

3.  To  conduct  a  useful  campaign  against  such  a  large  fire 
as  an  "experiment  of  opportunity"  calls  for  pre¬ 
planning:  an  appropriate  program  is  outlined  in  Part  I. 

One  issue  which  has  not  been  addressed  in  FY  1985  is  the 
difference  in  smoke  characteristics  between  forest/agricultural 
and  urban/industrial  fires,  i.e.,  between  wood  and  oil  or 
plastics  as  fuels.  While,  for  example,  in  the  SCOPE  briefing 
at  the  National  Academy  of  Sciences  on  13  September  1985  there 
is  repeated  reference  to  the  black,  sooty  smoke  from  urban/ 
industrial  fires  as  compared  with  the  white  smoke  from  burning 
wood  (forest  fires),  above  the  condensation  level  in  the  atmos¬ 
phere  there  may  be  little  difference:  certainly  the  oil  fire 
plume  of  Fig.  B-l  does  not  look  much  darker  than  the  neighboring 
water  or  ice  clouds.  It  may  be  possible  to  address  late-time 
soot  characteristics,  in  particular,  hygroscop ic i ty  and  effec¬ 
tiveness  as  condensation  nuclei  between  wood  and  oil  fires.  We 
know  that  f resh  soot  is  hygrophobic,  and  laboratory  experiments 
probably  cannot  study  the  properties  of  aged  soot.  However,  if 
one  can  find  satellite  imagery,  in  particular  mu  It ispect ra  1 
LANDSAT  (high-resolution)  imagery  of  oil  and  forest  fires,  this 
may  provide  quantitative  information  on  the  particle  growth  and 
change  in  optical  properties  over  times  of  hours  to  days.  We 
hope  to  begin  to  address  this  issue  in  FY  1986. 
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APPENDIX  C 


ATMOSPHERIC  MOISTURE  AND  ITS  VARIABILITY 

Figure  C-l  shows  the  atmospheric  water  vapor  mixing  ratio 
as  a  function  of  altitude,  combining  data  from  the  U.S.  Standard 
Atmosphere,  1976  with  some  "representative"  results  for  Munich 
in  January  and  July  1985  (kindly  provided  by  Dr.  J.K.  Angell, 
NOAA/ARL ) .  Figure  C-2  lists  the  relative  humidity  as  a  function 
of  altitude  for  the  U.S.  Standard  Atmosphere  (mid-latitude)  and 
for  tropical  and  subarctic  winter  models  from  the  IAMAP  radia¬ 
tion  atmosphere  (see  McClatchey  et  al.,  1979,  essentially  an 
AFGL  LOWTRAN  model). 

At  mid-latitudes,  the  relative  humidity  in  the  troposphere 
is  typically  50  percent  on  the  average;  in  a  subarctic  winter 
situation  the  relative  humidity  is  often  higher  in  the  lower 
troposphere  because  the  air  mass  has  cooled  but  not  lost  any 
moisture  by  precipitation.  In  the  tropics,  the  relative  humi¬ 
dity  is  low  above  the  boundary  layer  (where  moisture  is  picked 
up  from  the  ocean)  except  near  the  (very  cold)  tropical  tropo- 
pause  which  serves  as  a  "cold  trap"  for  limiting  the  amount  of 
moisture  that  can  enter  the  (very  dry)  stratosphere. 

Note  that  the  atmospheric  "cold  trap"  works  dynamically, 
not  statically.  Rising  thunder  clouds  are  warmed  by  condensing 
excess  moisture,  which  creates  buoyancy  that  enables  the  cloud 
to  rise  higher.  The  ultimate  height  limit  is  set  by  the  very 
low  ambient  tropical  tropopause  temperature  and  the  stability 
of  the  lower  stratosphere.  (For  an  experimental  investigation 
of  the  transport  of  water  vapor  by  tropical  thunderclouds  see 

IMaif  ah  a  1  flfl«91  .1 
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FIGURE  C-l 


Variability  of  atmospheric  water  vapor. 

(USSA  1976  and  Munich  January  and  July  1985) 
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FIGURE  C-2.  Relative  humidity  as  a  function  of  altitude  for 

IAMAP  radiation  atmospheres.  (McClatchey  et  al., 
1979. ) 
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Table  C-l  lists  precipitable  water  integrated  from  the 
surface  to  500  mb  (~  5.5  km)  at  selected  locations. 2  Note  that 
the  total  precipitable  water  varies  on  the  average  by  a  factor 
of  ten  between  different  locations  and  seasons:  at  high  lati¬ 
tudes  in  winter  (i.e.,  at  very  low  temperatures)  the  moisture 
content  is  very  low  (<  0.5  g/cm  ) ;  at  middle  latitudes  in 
summer,  the  moisture  is  moderately  high  (~  2  g/cm^)  while  in 
the  humid,  maritime  tropics  the  moisture  content  is  very  high 
( -  4  to  5  g/cm^ ) . 

Smaller  scale  structure  in  the  water  vapor  column  can  be 
obtained  from  various  satellite  sensors: 

(a)  McMurdie  and  Katsaros  (1985)  have  analyzed  the  inte¬ 
grated  water  vapor  column  obtained  from  the  scanning 
multichannel  microwave  radiometer  (SMMR)  on  SEASAT 

for  a  North  Pacific  cyclone  on  10  to  12  September  1978. 
They  find  that  the  water  vapor  column  varies  by  a 
factor  10  to  20  as  one  goes  through  the  surface  front 
of  this  weather  system.  Presumably,  the  same  effect 
is  observed  on  land,  but  a  sea  surface  provides  a 
particularly  benign  background  for  SMMR. 

(b)  Even  smaller  scale  structure  is  shown  on  the  "WV" 
channel  (5.7  to  7.1  ym)  which  is  present  on  the 
European  geostationary  meteorological  satellite 
METEOSAT  and  also  on  current  U.S.  GOES  geostationary 


^ At  mid-latitudes  thunderclouds  sometimes  rise  well  into  the 
stratosphere,  so  that  they  can  serve  to  transport  water  vapor 
to  those  altitudes.  It  is  not  known  how  effective  this  mech¬ 
anism  for  water  vapor  transport  is,  as  some  of  the  moisture 
will  go  back  into  the  troposphere  as  the  cloud  dissipates. 

^Standard  Rawinsonde  measurements  of  dew  point  (humidity) 
work  for  temperatures  above  -20  to  -30  C  (~  400  to  500  mb 
or  5  to  7  km).  At  lower  pressures  (higher  altitudes)  the 
air  is  so  dry  that  reliable  measurements  of  relative  humidity 
cannot  be  obtained  from  a  Rawinsonde,  and  in  any  case  the 
contribution  to  integrated  precipitable  water  is  small.  Note, 
however,  that  the  data  base  for  these  high  altitudes/low 
humidities  is  not  good. 
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500  mb)  AT  SELECTED  LOCATIONS: 


weather  satellites.  The  main  contribution  to  this 


channel  comes  from  a  layer  approximately  320  mb  thick 
which  lies  someplace  between  the  800  and  250  mb 
levels  (2  to  11  km),  depending  on  temperature  and 
humidity.  As  an  example  of  this.  Fig.  C-3  shows 
roughly  comparable  IR  (10.5  to  12.5  um)  and  WV 
(5.7  to  7.1  um)  views  of  the  earth.  The  IR  channel 
shows  the  thermal  emission  from  clouds,  while  the  WV 
channel  also  gives  a  return  from  water  vapor  in  the 
atmosphere.  Both  signals  show  64  gray-shades,  i.e., 
the  ratio  of  maximum  to  minimum  intensity  that  can  be 
obtained  from  digital  imagery  can  be  as  much  as  64, 
although  the  photographic  image  does  not  have  this 
wide  a  dynamic  range. 

The  conclusion  is  that  atmospheric  moisture  is  highly 
variable,  as  a  function  of  the  "origin"  of  a  particular  air 
mass  between  the  extremes  of  a  cold,  dry  continental  polar  (cP) 
air  mass  and  a  warm,  wet  maritime  tropical  (mT)  air  mass.  The 
horizontal  scale  of  a  midlatitude  air  mass  is  of  the  order  of 


1000  km,  so  that  there  can  be  a  significant  variation  in  atmo¬ 
spheric  moisture  at  different  locations  on  a  particular  day, 
perhaps  by  an  order  of  magnitude. 


APPENDIX  D 


EFFECTS  OF  CLOUD  NONUNIFORMITY 


When  a  given  amount  of  smoke  is  injected  into  the  atmo¬ 
sphere,  the  cloud  will  hardly  ever  be  distributed  uniformly  in 
a  horizontal  plane  (see  Fig.  2).  The  effect  of  this  will  be 
that  with  a  non-uniform  distribution  (Case  2)  more  radiation 
will  always  be  transmitted  through  the  cloud  than  for  a  uniform 
distribution  (Case  1),  just  because  of  the  nature  of  the  (expo¬ 
nential)  extinction. 

Figure  D-l  illustrates  this  for  one  example.  The  total 
amount  of  smoke  corresponds  to  an  extinction  optical  thickness 
X  when  it  is  distributed  uni  f  ormly--in  Case  1.  In  Case  2,  the 
same  total  amount  of  smoke  is  divided  between  two  regions:  a 
region  of  area  Athin  in  which  the  optical  thickness  is  smaller 
by  a  factor  F,  i.e.,  is  X^/F,  and  another  region  of  area  A  t  ^  i  c  k 
in  which  the  optical  thickness  has  the  value  F  X^,  i.e.,  is 

comparably  greater.  The  relative  areas  Aj-^ick  ar>d  Athin  are 
determined  by  the  condition  that  the  total  amount  of  smoke  is 
the  same  as  in  Case  1,  i.e., 


Athick  ^  *1  +  ^thin  ^i/F  ~  (A^hick  +  ^thin^  ^1  (D.l) 

Figure  D-l  shows  the  effective  transmittance  ratio,  ^2^1 

through  the  cloud  as  a  function  of  optical  thickness  X ^  for 
representative  values  of  F  from  1  to  16,  where  the  ratio  At-^ick- 
Athin  is  Qiven  by  Eq.  (D.l).  Evidently,  for  F  =  1  there  is 
no  non-uniformity  so  that  T2/T^  =  1,  while  as  both  X^  and  F 
increase,  the  ratio  T2/T1  increases.  Thus,  e.g.,  for  X^  =  3, 

the  ratio  T2/T1  =  1.52  if  F  =  1.41  or  7.56  if  F  =  4. 


MTIO  AT  OPTICAL  THICKNESS  OF  TWO  KE6 tONS  IN  CASE  2 


